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ON  OBSERVER  PROBLEMS  FOR  SYSTEMS 
GOVERNED  BY  PARTIAL  DIFFERENTIAL  EQUATIONS  * 


J.S.  Barast  and  A.  Bensoussanft 

tUnivcrsity  of  Maryland 

tt University  Paris  Dauphine 
and 
INRIA 


1.  Introduction 


The  problem  of  observers  has  been  introduced  in  the  control  litera¬ 
ture  by  D.  Luenberger '|15LV  Let  us  consider  a  dynamic  system  which  is 
deterministic,  but  whose  inital  state  is  unknown.  An  observer  is  a  model 
which  mimics  the  behavior  of  the  physical  system,  and  in  particular  its 
state  becomes  closer  and  closer  as  time  evolves  to  the  state  of  the  physi¬ 
cal  system  There  is  a  great  deal  of  freedom  in  such  a  design  and  it  is 
important  to  investigate  various  kinds  of  observers. 

Since,  after  all,  the  observer  prblem  presents  analogies  with  the  filter¬ 
ing  problem  (estimating  the  state  of  a  stochastic  dynamic  system),  although 
there  are  no  stochastic  disturbances,  it  is  natural  to  exploit  the  analogy. 
This  idea  has  been  used  by  J.S.  Baras  and  P.S.  Krishnaprasad  (4)  and  leads 
to  an  observer  which  is  different  from  Luenberger's  observer.  It  presents 
several  advantages.  In  particular,  it  is  obtained  in  a  constructive  way  and  it 
has  robustness  properties. 

From  the  very  definition  it  applies  identically  when  there  are  distur¬ 
bances,  whereas  the  Luenberger  observer  is  strictly  limited  to  the  determin¬ 
istic  case  and  is  not  obtained  in  a  constructive  way.  Also  it  may  apply  to 
more  general  cases,  in  the  sense  that  when  the  Luenberger  observer  exists, 
the  observer  arising  from  the  Kalman  filter  theory  exists  as  well. 

In  this  article  we  consider  dynamic  systems  whose  evolution  is  governed 
by  a  parabolic  partial  differential  equation,  or  more  generally  a  differential 
operational  equation  in  the  sense  of  J.L.  Lions  [121 

The  Luenberger  theory  has  been  extended  to  infinite  dimensional  sys¬ 
tems  (see  in  particular  M.J.  Chapman  and  A.J.  Pritchard  [61,  A.  Ichikawa 
and  A.J.  Pritchard  ( 10)).  We  explore  here  the  observer  based  upon  Kalman 
filter  theory. 

Let  us  discuss  here  an  example  to  present  the  main  reults  of  the  paper. 

Let  //  be  a  smooth  domain  of  R".  Consider  the  P.D.E. 
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S 


dK  A 

a?"A> 


0  in  /? , 

0. 

>(*.  o)  «  y0(x)  e  H  -  L2(/2)  , 


a^lr 


(1.1) 


where  F  =  df?  denotes  the  boundary  of  17  The  state  ye  is  unknown.  If  in¬ 
stead  of  (1.1)  we  were  to  have 


-  A y  +  Ay  -  0  ,  A  >  0  , 


& 
dt'  i 


0. 


(1.2) 


y(*>  o) 


then  the  system  is  stable  and  y(x,  T)  ->  0  as  T  -♦ <».  In  other  words,  whatever 
ya  is,  the  state  comes  closer  and  closer  to  0  as  T  evolves  and  thus  becomes 
“more  and  more  known.”  An  observer  could  be  simply  the  model  (1.2)  itself 
with  an  arbitrary  value  of  y0  independent  of  the  available  observation. 

Such  a  stability  property  is  not  present  in  the  model  (1.1).  Let  us 
assume  that  we  observe 


y\. 


(i.i)' 


i.e.,  the  value  of  the  state  on  the  boundary. 

An  observer  in  the  spirit  of  Luenberger  would  be  the  following  model 


dm 
d  t 


-  Am  -  0  in  Ft  , 

dm  .  . 

57  -  (*  -  "Or  . 


(1-3) 


m(x,  o)  -  mc(x)  , 

where  ma  is  arbitrary.  The  error  rj  «  y  -  m  appears  as  the  solution  of 


~  A7  -  0  in  Q  , 

I 

-  0, 


dri 

di  +  i 


(14) 


TJ{x,  o )  -  7]0{x)  -  y0  -  ma  . 
Multiplying  (4)  by  T]  integrating  over  O  yields 

l  d 

r 

a 
But 

J  \Drj\ 2  dx  +  Jr?  dr  Zfijrfdx,  0  >  0  . 


5"  2i  J  +  J  dx  +  J  n7  dr  -  0 

a  a  r 


n 


a 


iii 
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and  thus 


J  TJ7(x,  t)dx  <  Ce~2fil  , 
n 


which  proves  the  exponential  decay  of  the  error. 

The  theory  developed  in  this  paper  leads  to  the  following  observer: 


jj  -  Am  =  J  P(x,  £  r)(r(£,  r)  -  m(£  t))  d£  , 


§5i-- «• 


(1.5) 


— — ^  m(x>  ° *  “  ’ 

where /?(*,  £,  /)  appears  as  the  solution  of  a  Riccati  equation,  connected 
to  a  filtering  problem,  or  by  duality  to  a  control  problem.  We  study  the  type 
of  control  problems  which  may  be  introduced  in  order  to  derive  exponen¬ 
tial  decay  for  the  error. 


2.  Setting  of  the  Problem 
2 . 1  Notation  and  Assumptions 

Let  V,  H  be  two  separable  real  Hilbert  spaces,  such  that  identify  ing  H 
and  its  dual  H\  one  has 


V  c  H  «  H'  c  V  (2.1) 

each  space  being  dense  in  the  next  one  with  continuous  injection. 

We  denote  by  ((  ,  ))  II  II  and  (  ,  )  |  |  the  scalar  product  and  norm  in 
V',  H,  respectively.  For  another  Hilbert  space  X  we  shall  use  the  notation 
(  ,  )x  and  |  |x .  We  denote  by  <  ,  >  the  duality  V,  V’.t 
Let  /4(r)  be  a  family  of  operators  such  that 

A(  )  e  Lm(o,  V;  V')) 

(2.2) 

<A(r)0.  0>  +  A|0i2  £  or H2  ,  V  0  €  V  , 

a  >  0,  /iSO.  V  t . 

We  consider  a  dynamic  system  whose  evolution  is  governed  by  the  dif¬ 
ferential  operational  equations  (see  J.L.  Lions  (12)) 


%  *  MOy  -  / . 
yio)  -  ya  . 


(2.3) 


t  More  generally,  <  ,  >*  will  represent  the  duality  between  a  Hilbert  space 
X  ands  its  dual  X’  and  Ax  the  canonical  isomorphism  between  X  and  X’. 


I 
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where 


/  e  Li(o,  oo;  V")  given  , 
ya  e  H  ,  yB  unknown  . 


(2.4) 

(2.5) 


It  is  well  known  that  for  any  y0,  (2.3)  defines  the  state  >•(•)  in  the  sense 


2  “•  ,A  -  e  ljoc(o,  «;  VO 


>  e  Hoc(°.  °°;  vo .  dt 
We  perform  an  observation  on  the  state  y(  )  as  follows: 

2(0  -  C(0y(0. 

where 

C(  )  e  L*°  (o,  «>;  5C(V;  F))  , 

F  being  a  given  Hilbert  space. 


(2.6) 


(2.7) 


2.2  Definition  of  an  Observer 

An  observer  is  a  function  m(t)  measurable  with  values  in  H,  whose 
value  can  be  computed  at  each  time  t  in  terms  of  the  known  data  (in  parti¬ 
cular  the  observation  *(■))  and  such  that  the  error 

e(t)  -=  y(t)  -  m{t) 


satisfies 

KOI  -♦  0  as  t  -*  «> . 

Therefore  m(t)  will  reasonably  estimate  the  state  of  the  system  at  time 
t  (at  t  -  0  this  estimate  may  be  very  bad,  but  it  improves  more  and  more  as 
r  -» oo ;  it  is  of  course  nice  to  get  an  exponential  decay  for  |e(OI). 

This  problem  has  been  extensively  studied  in  the  finite  dimensional 
case,  starting  with  the  seminal  work  of  D.  Luenberger  115).  For  infinite 
dimensional  systems,  the  Luenberger  observer  has  been  exended  in  a  natural 
fashion,  however  the  research  has  mostly  concentrated  on  the  design  of  com¬ 
pensation  (i.e.,  use  the  possibility  of  controlling  the  system,  so  that  the 
global  system  made  of  the  system  itself  and  of  the  observer  is  stable). t 
While  the  control  is  generally  present  in  such  contexts,  it  is  worthwhile  to 
separate  the  observer  problem  from  the  compensation  problem,  and  consi¬ 
der  cases  similar  to  (2.3),  (2.7)  when  there  is  no  control.  In  particular, 
nothing  can  be  done  to  stabilize  the  system  itself. 

We  shall  thus  prefer  to  introduce  observers  constructed  in  a  different 
way,  following  ideas  introduced  by  J.S.  Baras  and  P.S.  Krishnaprasad  [4], 
which  consist  in  artificially  randomizing  the  problem  and  using  Kalman  fit¬ 
ter  theory  for  infinite  dimensional  systems. 


t  For  details  in  this  direction  see  MJ.  Chapman  and  A.J.  Pritchard  16], 
A.  Ichikawa  and  AJ.  Pritchard  (10). 
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2.3  /4  Randomizing  System 

The  theory  of  linear  stochastic  infinite  dimensional  systems  has  attracted 
a  considerable  interest  in  the  literature.  One  of  the  main  objectives  has  been 
to  develop  a  rigorous  theory  of  the  Kalman  filter  applicable  to  distributed 
parameter  systems.  Among  the  main  contributions  one  can  refer  to  A.V. 
Balakrishnan  (1],  [2],  A.  Bensoussan  [5],  R.F.  Curtain  |8],  R.F.  Curtain 
and  A.J.  Pritchard  |9]. 

Let  us  describe  briefly  what  we  shall  need  for  our  purpose.  Let  (/?,  3,  P) 
be  a  probability  space  equipped  with  a  filtration  satisfying  the  usual  con¬ 
ditions. 

Let  £  be  a  Hilbert  space,  a  generalized  E  Wiener  space  is  a  stochastic 
process  indexed  by  an  element  e.(  )  €  l^oc(o,  £'),  denoted  by  ft,"  (.oS) 
satisfying 


/rf*  is  a  7'  Wiener,  V  e.  , 

(2.8) 

r i 

em'*p‘,'-I  <*•&).  a~e  *•&)><&, 

o 

(2.9) 

the  map  e.  -»//**  is  linear  . 

(2.10) 

We  assume  that  such  a  generalized  E  Wiener  process  exists.  Similarly 
we  shall  assume  the  existence  of  a  generalized  F  Wiener  process  vh,’  and  we 
suppose  that 

//'*  and  v,  are  independent  .  (2.11) 

Let  C(  )  e  Z~(o.  «;  V'))  and  M{)  €  C (o,  «»;  5 6(£;  £)),  we  shall 

consider  the  stochastic  proceses  ft,  ”  and  \\  *  indexed  respectively  by 
v(  )  6  L,Joc(o,  »;  VO  and  /.(•)  €  tfoc{o,  «o;  F').  They  are  clearly  indepen¬ 
dent  ?'  processes  and  one  has 

•  •  I A  S 

EfP,  y,fl,  ’*  -  f  <GA~e'G*vx(A),  v2(/l)>  dX  .  (2.12) 

O 

Finally,  let  there  be  a  family  of  random  variables  indexed  by  h  e  H  denoted 
by  such  that 

V  h  ,  is  a  Gaussian  with  mean  O,  and 

v  (213)  , 

-  iPeTh,h),  P„  €  KH;H).  >v 

symmetric,  semipositive  definite. 

is  independent  from  fi'  and  »'{* .  (2.14) 

Let  r\t ,  j),  r  £  i  be  the  Green's  operator  corresponding  to  A(r),  i.e., 

rxt,  s)  €  Z(H,H)  and  a(r)  -  F\t.  s)h  (2.15) 

is  the  solution  of 


da 

dt 


+  Aa 


0  ,  t  >  j,  a(j)  -  . 
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For  simplicity  we  write  J~[t)  -  r\t,  o). 

A  stochastic  linear  system  is  a  process  y,  indexed  by  h  e  H,  defined  by 
the  formula 

>•?  -  (y(0,A 0  +  £r’(0*  +  /'?VV)*  ,  (2.16) 

where  y(r)  is  the  solution  of 
d\ 


J7  +  ^Oy  *  /(0  .  y(°)  -  y0 

y0  given  in  H  . 


(2.17) 


We  are  to  be  careful  in  the  interpretation  of  //f  r  in  fact, 

consider  5  ->  G*(s)£*(r.  j)/i  for  s  <  r;  it  coincides  with  G*(s)/?(s),  where  /? 
is  the  solution  of 

ff  +  A*(s)fi  «  0  ,  AO  -  A  , 

and  AJ)  €  L?{o,t,\r),  hence  s  -» G*(s)T*(t,s)h  «=  e'.fs)  belongs  to  L2(o,r;  £’). 
Since//**  depends  only  on  the  restriction  of  e.  on  ( o ,  t),  the  quantity//**' 
is  well  defined  for  any  /. 

We  define  the  observation  at  time  T  as  follows.  Let  /.  e  Lfoc(o/,  7££') 
and  71  be  the  solution  of 


-fi  +  A'n  «  c7.  , 

*7(0  -  o. 


Then  one  sets 
.7 


zfj  -  J  </.(/),  C(r)y(0>  +  £’?(c)  +  //?  7  +  ^  A 


(2.18) 


(2.19) 


£ 


2.4  Kalman  Filler 

This  problem  can  be  stated  as  follows.  Let 

ZT  «  <7(/\/.  €  mo,  T,  £•))  . 

Find 

y^  -  £[y*  I  Z  J  .  (2.20) 

Without  redoing  the  theory,  for  which  we  refer  to  A.  Bensoussan  |5J,  we 
shall  only  recall  that  it  can  be  obtained  by  solving  a  deterministic  control 
problem  related  to  maximum  likelihood.  We  assume  that  M  is  invertible  as 
well  as  Pg  and  set 

R  -  €  *<£,£').  (2.21) 

We  introduce  the  control  problem 


y(°)  -  y0  +  £  • 


in  which  £  and  e(  )  arc  «he  decision  variables. 

We  are  interested  in  minimizing  the  following  cost: 

7(£,  e( ))  -  /  (</i£e,  e>  +  <R(£ -  Cy),  f-  C>>]  <* 


(2.23) 


where  on  the  right  side  of  (2.23)  the  function  £(■)  is  given  in  Lj(o,  T;  F). 
Define 

D,  -  GA~E  G*  , 

(2.24) 

D2  -  C  RC  , 

and  consider  the  pair  y(),  p( •)  defined  by  the  system  of  coupled  equations 

^  +  A(t)y  +  D,(0/>  =  /. 

-  ^  +  >4 *(/)/?  -  Z)j(r)y  -  -C*R£ , 

>(»)  -  y*  -  **«,/>(<>)  . 
p{T)  =  0. 

Then  the  optimal  control  is  given  by 

e(r)  =  -v^Vw/*)  , 

£  -  -P0p{o)  . 

The  decoupling  theory  leads  to  the  Riccati  equation  (written  formally) 

4  AP  4  PA*  4  PD2P  -  D,  , 


(2.25) 


(2.26) 


(2.27) 


P(o)  «  P0  . 


and  the  linear  equation 


^4*r  -  /  4  PC*R{£-Cr), 
>io)  -  yc  . 

It  can  be  proven  that  if  e  Li(o,  T ;  F")  is  defined  by 

T 


(2.28) 


J  </..C>dt  -  (r(T)  -  y(T),  h)  ,  V  f  , 


then  the  quantity  (2.20)  is  given  by 


>*  -  /  *  +  O^).  A) . 


3.  Study  of  the  Operator  P(T) 
3.1  Definition  of  P(T) 

Consider  the  coupled  system 


+  A6  +  D\fi  -  0  , 
A*fi  -  D76  =  0  , 


<*(0)  -  -/’«#*). 


Then  we  set 


An  -  *• 


P(T)h  =  -d(T) 


The  system  (3.1)  is  related  to  the  following  control  problem: 

+  A{t)a  =  Ce  , 

«(<?)  -  (3.3) 

(#.«*(•))  -  +  J  lklJ  +  <D2a,o'>]  </;}  +  (A,  a(J)) . 

O 

We  deduce  easily  that 

T  (nT)h.  h)  -  -  inf  ./?(£  «(•))  •  (3.4) 

*.«(•> 

We  can  also  characterize  this  quantity  in  a  different  way.  Consider  the  con¬ 
trol  problem 


A*fi  -  CV.  .  *.  e  L2(o,  T;  F')  , 

PW)  - 


and  the  cost 


#(*.(•»  -  (P0P{o).fco))  +  J  +  J  <Dl0,p>dt. 


Then  one  has  also 


(P(T)h,h)  -  inf /£(*.(  ))  . 


1P1 


-9- 


3.2  Detectability  and  Stabilizability 

definition  3.1.  We  shall  say  that  the  pair  A(),  C(  )  is  detectable  if  V  h, 
3  <pT.'h  such  that 

A  +  /V’*l l  dt  <  Kh  (3.7) 

O  O 

T  h 

independent  of  T,  where  ft  '  represents  the  solution  of  (3.5)  corresponding 
to  <p\-h.  □ 

In  the  stationary  case,  i.e..  A,  C  independent  of  time,  it  is  sufficient 
to  assume  the  following 

definition  3.2.  We  say  A *,  C*  is  stabilizable  if  V  h,  3  e  L?(o,  «■;  F') 
such  that  the  solution  y*  of 

dyh  a*  j, 

-f—  +  A  yh  =  C  0.  , 


y*'{p)  =  h  , 


satisfies  yh  e  L2(o,  ~>;  H).  □ 


proposition  3.1.  In  the  stationary  case,  if  A*,  C*  is  stabilizable,  then  A, 
C  is  detectable. 

Proof.  Define 

*•’  (0  -  4>.{T-t),  p-\t)  -  ^(F,  t)  . 

T  h  T  h 

Then  clearly  /?  '  (/)  is  the  solution  of  (2.5)  corresponding  to  <p. '  (r),  and 

j  i*r*i \\.dt  -  j  \£\ 7  dt  <  j  wh.\ 1  di, 

o  o  o 

JT]/f‘h\7dt  -  /Vi*  dt  <  f  Ir*i 2  dr, 

o  o  o 

and  the  desired  property  follows.  D 

theorem  3.1.  Assume  that  A(  ),  C(  )  is  detectable,  then 

I WIm //,//)  -  P  ■  (3-8) 

Proof  It  follows  from  (3.5)  that 

T  T 

r  W<>)\2  4  J  <A  V.  P>  dt  -  j-  |Ap  +  J  </?,  CV.  >  dt  . 


Hence 


j-  po)\7  4  of  dr  £  j-  |/ip  4  J  <P  C>.  >dt  4  Aj  ]/}•}  dt 


and  ihe  detectability  condition  implies  immediately 

)f'h(o)\  <  k;  ,  jVT  V  <  K' . 

o 

Therefore  from  the  definition  (3.6)  one  has 

(P(T)h,  h)  <  ./*«>•'*)  <  K-- . 

This  and  the  fact  that  P(T)  is  symmetric  positive  semidefinite  implies  the 
desired  result.  D 

theorem  3.2.  If  there  exists  a  family  7~{  )  e  Lr  (o,  <*;  £(5;  l"))  such  that 

<(A(i)  +  mC(t)l  <?>  £  (3.9) 

then  the  pair  A(),  C(  )  is  detectable. 

Proof  Indeed  pick  the  feedback 

<t>.  -  -r*p. 

The  corresponding  trajectory  is  given  by 

-  ^  4  (A*  4  C*r*)p  =  0  .  P(T)  =  h  , 

and  thus  Jj  \\p\?di  is  bounded  by  a  constant  independent  of  T.  The  de¬ 
sired  result  follows.  □ 

3.3  Invertibility 

We  turn  now  to  the  question  of  invertibility  of  the  operator  P(T).  We 
shall  need  another  property.  Consider  the  dynamic  system 

da  .  ~ 

4  Aa  *  Ge  , 

a  t 

(3.10) 

a(o)  - 

definition  3.3.  We  shall  say  that  the  pair  A(),  G()  is  controllable  if 
V  h,  3  eT  h  and  such  that 

\f*t  4  /  |eT-*J2  di  4  /  \aT  h\2  di  <  Lh  (3.11) 

•  0 

and 

aT  h(T)  -  h  ,  (3  12) 

where  aT •*  designates  the  solution  of  (3.10)  corresponding  !;  ’h  and  e1  •*. 
The  constant  Lh  is  independent  of  T.  D 
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We  can  give  an  example  of  controllability  (probably  the  only  one  really 
implementable). 


proposition  3.2.  Assume  that 

D\  is  invertible  . 

Then  the  pair  A(-),  G( )  is  controllable. 

Proof  Let 


(3.13) 


d<*> 


<Ph  F  L<(o,  °°.  10  .  -Jf  e  Li (o,  «■;  V')  ,  <?h{o)  =  h 


We  set 


aT  h(‘)  =  <Ph(T-  t),  ?■*  =  aT.>(o)  =  ph(T) 


T.h 


and 


«T>*(0  =  /lf'c^DKO'1  +  A(r)ffT»(,)] 

Then 

£T’V  «  Vh(T)f  =  2jr(^,^)  <*  +  | h? 

O 

<  w’  *  2/ "{«>,.  jr)  <*. 

O 

rT  rT  r T 

J  \aT  h?dt  =  J  k?/,(T-r)!2  dt  =  j  k?,(r)|J  A  <  J  | <?h?  dt  . 

0  o  00 

and 

/V-*(0P*  <  c[jr||^||2vyr  +  /TBflr7-.*(0H2* 

<  C 

Finally  from  (3.4)  we  obtain 


(3.14) 


f  llsrllv*  */"«<». *’<* 


.  w  /yT ,  /1 

G(t)eT-h(t)  =  -777 —  +  A(/)flrr-*(/) 


and 


ar'*(T)  -  /j  . 


The  desired  result  has  been  established.  0 


remark  1.  As  is  well  known  (cf.  R.  Lattes  and  J.L.  Lions  (11)),  we  cannot 
solve  a  priori  the  backward  problem 


da 


j  +  Aa  -  Ge  ,  a(T)  -  /i  . 


(3.15) 


The  situation  is  different  from  ordinary  differential  equations.  The  pro¬ 
blem  (315)  (for  given  e)  is  a  priori  ill  posed.  We  refer  to  J.L.  Lions  [14] 
for  a  detailed  study  of  these  ill  posed  problems  in  the  context  of  control 
theory.  □ 

We  shall  now  consider  control  problems  similar  to  those  described  in 
Section  3.1.  The  response  is  described  by 


da  .  ~ 

3—  4  Aa  =  Gc  , 
dt 

a(o)  =  g  . 


(3.16) 


We  impose  the  constraint 


a(T)  =  h 


(3.17) 


and  minimize  the  cost  (recall  that  Pa  is  invertible  (cf.  Section  2.4)) 

T 

*}(£«(•))  -  (/>;'££)  +  J  \\e\l  +  <D2a,  a>]  dt .  (3.18) 

O 

Note  that  (3.17)  must  be  considered  as  a  constraint  and  not  an  initial 
condition. 

theorem  3.3.  Assume  that  the  pair  A( •),  G()  is  controllable.  Then  the  con 
trol  problem  (3.16),  (3.17),  (3.18)  has  a  unique  solution.  There  exists  a 
unique  pair  5,  /?  such  that 

—  +  Aa  4  D\fi  =  0  , 

-  ~p-  4  A  P  -  Dfa  =  0. 


B(o)  ■=  -Pj{o)  , 
5(T)  -  h  , 


(3.19) 


and  the  optimal  control  is 


e(z)  -  -AE'G\t)p(t)  . 
I  -  -Pofi(o)  . 


(3.20) 


Proof  We  follow  the  technique  introduced  by  J.L.  Lions  [13]  to  deal  with 

this  type  of  ill  posed  problem.  We  shall  penalize  the  constraint  (3.17). 

^  -  _  .  _  _ 

a.  Existence  and  Uniqueness  of  the  Optimal  Control 

Let  rh  -  {e(),  £\a{T)  -  A}.  By  the  controllability  assumption  Fh  is 
not  empty.  It  is  a  convex  closed  subset  of  the  Hilbert  space  Li(o,T;  E)  x  H 
and  the  functional  (3.18)  is  a  coercive  quadratic  functional.  Hence  the  exis¬ 
tence  and  uniqueness  of  the  optimal  control  ?(•),  £. 


b.  Necessary  and  Sufficient  Conditions  of  Optimality 


Consider  P0  (i.e.,  rh  with  h  »  0).  It  is  obviously  a  closed  subvector 
space  of  L2(o,T;  £)  x  H.  Noting  that  ?(•)  -*•  Ae(),  £  +  A4  V  *(■),  g  be¬ 
longing  to  F0,  we  deduce  easily  from  the  relation 

4(1  +  ^.  )  +  M  *  4(£e(» 

that  the  following  condition  must  hold 

J/  r  T 

£#)  +  J  <AEe,e>dt  +  J  <D275,a>dt  =  0,  (3.21) 

0  0 


and 


V  £,  e(  )  in  £c  , 

a  being  the  corresponding  solution  of  (3.16), 
75  designates  the  optimal  trajectory  . 


This  condition  is  sufficient,  since  V  £,  e(  )  e  £A,  £  -  e{  )  -  e(  )  e  £e. 
Hence 

✓ 

7  T 


<(f.  «(•»  -  (p 


and 


J  </1£e,  e>  dt  +  J  <D2a.  a>  dt 
0 

*())  >  #f(£  ?(•))  ■ 


c.  Adjoint  System.  Uniqueness 

Assume  that  we  have  a  solution  of  (3.19).  Then  define  e( •)  and  .5  by 
(3.20).  Let  us  show  that  it  is  optimal.  Clearly  75  is  the  trajectory  corre¬ 
sponding  to  the  control  ?(■),  £  and  it  is  admissible  (belongs  to  rh).  It  is 
easy  to  check  that  it  satisfies  (3.21).  It  is  thus  the  optimal  control  and 
therefore  ?(•).  £  and  75{)  are  uniquely  defined.  It  is  also  the  case  of  (i\o) 
and  -(dfl/dt)  +  A* fi.  This  implies  that  is  unique  as  well. 

d.  Adjoint  System.  Existence 

It  remains  to  prove  the  existence  of  the  pair  75,  P  solution  of  (3.19) 
The  penalty  technique  is  now  used.  Consider  the  functional 

^•f(f,e(  ))  -  #£(£«(•))  +  f  \a{T)  -  h?  .  (3  22) 

Then  the  control  problem  (3.16),  (3.22)  becomes  classical,  and  there  exists 
a  pair  ac,  pc  solution  of 


and  the  optimal  control  is 


ec  -  ~a~e'g*Pc  .  4c  -  -p0PA°)  ■ 

Now  from 

4'V.*(-))  ^  Kkj(f  \eTI')  <  C, 


we  deduce 


rT  rT 

J  K\2e  dt  J  <D2ff£,  0'f>  dt  <  C* 

F  \ac(T)  -  h?  <  Ch  , 


(3.24) 


and  since  ac  is  the  trajectory  corresponding  to  ec  and  £c  we  have 

/T  H5riir<i'  *  ^  <  c>.t 


(3.25) 


(this  last  constant  may  depend  on  T).  Now  for  any  £,  e(  )  we  have  (neces¬ 
sary  condition  of  optimality  for  the  problem  (3.16),  (3.22)), 


(p0'4c’4)  +  f  <A£ec,e>dt  +  J  <D2ac,  a>  dt 


(3.26) 


♦  f  «*AT)  -  cdj))  -  o  . 


Let  us  pick  £  «  £T’*,  e  -  eT-*  where  *  is  arbitrary  in  H.  We  deduce 

3.24),  (3.25),  that 

•  *)  I  5  C*.7\*  *  C*  • 


from  (3.26)  and  the  estimates  (3.24),  (3.25),  that 
r<*e( T)  -  h 


since  in  this  context  h,  T  are  fixed.  From  the  Banach-Steinhaus  theorem  it 
follows  that  \(ac(T)  ~  h)le\  S  Ch  T.  Therefore,  also  using  (3.23),  one  has 

f&V  4  JT,A"'  *  <  C>.T  ■  (3.27) 

o  o 

With  the  estimates  (3.24),  (3.25),  (3.27)  we  can  extract  a  subsequence  ac, 
fic,  converging  weakly  to  (dWdi),  (dflldt),  in  Ljio,  T;  V').  We  obtain  easily 
that  H,  0  is  a  solution  of  (3.19),  the  proof  has  been  completed.  D 
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- .  1 < *y : w ' H*W> M •* -* *?i‘I K< 


theorem  3.4.  Under  the  assumptions  of  Theorem  3.3  one  has 

P(J)  -  -Q(T)h,  (3.28) 

where  Q(T)  e  %j(H;  H)  symmetric  positive  semidefinite.  Morever, 

\Q(T)\ t{H.H)  <  C.  (3.29) 

Proof  Consider  (3.23).  the  quantities  ac,flc  are  linear  continuous  functions 
of  h.  Hence  we  can  write 

pc(T)  -  -Qq)h  ,  (3.30) 

where  it  is  easily  seen  that  QC(T )  e  SC<//;  H)  symmetric,  positive  semide¬ 
finite.  Moreover,  one  easily  computes 


inf  4'c(4,  *(•))  -  (Qc(T)h,  h)  . 
i.'O 


(3.31) 


Since  the  left  hand  side  increases  as  £  decreases,  it  follows  that 
QC(T)  is  an  increasing  family  of  symmetric  positive  semidefinite  operators 
in  H). 


But 


Therefore 


4'r(£  e( ))  -  4(4,  *(•))  V  4,  e(  )  in  rh  . 


(Qc(T)h.h)  <  inf  4(<f,e(» 

#•«(■)«  rk 


<  4(4Tk,c?>')  <  c* 


(3.32) 


independent  of  T  and  c.  Necessarily  Qe(T)  T  Q(T)  e  'tlH,  H )  positive 
semidefinite  and 

\Q(T)\ *<„,„)  S  C  (3.33) 


independent  of  T.  But 

inf  4'Cl£>  *(  ))  "»  >nf  4(1  «0>  •  (3-34) 

£.*(■)  £.«<>«  rh 


Indeed, 


(K'tr +  ♦  f  <D2ac,a(>d,  +  '-\at(T)-h? 

£  (P’o'H)  *  /T|F|2  4fr  4-  f  dr  . 


and  from  the  weak  convergence 


•/*.  , 
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This  implies  the  strong  convergence  of  £c,  ec ,  ac,  to  £,  e,  ZF,  and  also  that 
( ]/e)\ac(T )  -  h|2  -»0.  t  Therefore  (3.34)  is  demonstrated  and  thus 

inf  /£(£*(•))  -  ( Q(T)h,h ).  (3.35) 

Moreover,  since  fic(T)  -*  fi(T)  in  H  weakly,  we  deduce  the  property  (3.28), 
too.  □ 

rim  ark  3.2.  Before  it  was  used  in  treating  ill  posed  problems,  the  penalty 
technique  (which  goes  back  to  Courant  (7)),  has  been  widely  used  as  an 
approximation  technique  for  the  control  of  systems  governed  by  partial 
differential  equations,  see  A.V.  Balakrishnan  [3],  J.L.  Lions  (13].  D 

We  can  now  compare  Q(T)  and  P(T). 

theorem  3.5.  Assume  that  A(),  C( )  is  detectable  and  A( •),  G( )  is  control- 
able  Then  one  has 


Q(T)P(T)  «  P(T)Q{T)  -  /.  (3.36) 

\Qn(H'H)  <  q,  \P(T)\tiH.H)  *  P  •  (3.37) 

where  p,  q  are  constants  independent  of  T. 

Proof  The  property  (3.37)  has  already  been  proven.  The  property  (3.36) 
follows  by  comparing  (3.1),  (3.2)  to  (3.19),  (3.28).  Indeed  in  (3.19)  set 
h  =  -P(T)h  Then  we  have  5  -  A  and  (by  uniqueness)  and 

JkT)  «  h  «  -Q(T)(-P(T)h)  -  Q{T)P(T)h  . 

A  similar  proof  is  made  to  prove  that  PQ  -  /.  D 

4.  Observer  Based  Upon  Kalman  Filter 
4.1  The  Model 

Motivated  by  the  form  of  the  Kalman  filter  see  (3.28)),  we  shall  de¬ 
fine  the  observer  by  the  solution  of  the  equation 

+  Am  -  /  +  PC*R(z-Cm),  (4.1) 

m(o)  -  mc 

where  mB  is  arbitrary  and  i  is  the  observation  corresponding  to  the  state  (2.3). 

The  writing  (4.1)  is  somewhat  formal,  since  P  is  not  defined  on  V'.  In 
fact,  we  shall  give  a  meaning  to 


y-m  -  rj  , 


(4.2) 


t  This  follows  also  directly  from  the  fact  that  ^a([T)-h)!c\  is  bounded  in  e 
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where  77  appears  as  the  solution  of 


j}  +  (A  *  PD2)tj  -  0  , 
*7(o)  -  ya  -  mo  ■ 


(4.3) 


The  solution  77(f)  of  (4.3)  is  defined  by  duality.  Indeed,  considering  (3.1) 
we  see  that  the  equation 


d/f 

~dT  *  (/*  +  DlP)P 


f(T)  -  h  , 


(4.4) 


has  a  solution  in  the  functional  space 

W>(0.  T)  =  {«>  €  L\oJ ■  V),$e  L2(o,T\  l'*),  /><9  e  Z.2(o,r;  .  (4.5) 


Then  the  value  tj(T)  is  defined  by 

W(T),h)  «  {/(o),y0-m0)  ,  V  h  ,  (4.6) 

which  defines  rj(T)  uniquely  in  H. 

Our  problem  amounts  to  studying  the  behavior  of  77(7")  as  T  -»  0 


4.2  Estimates 

We  shall  prove  the  following  (main)  result. 
theorem  4.1 .  Assume  A(  ),  C( )  detectable,  and 

<Di(f)v,  v>  2:  k\v\jf  V  v  e  V  .  t  (4.7) 

Then  one  has 

|y(r)  -  m(f)|  <  Cly0  -  mc|e'r'  ,  y>  0  .  (4  8) 

Proof  Considering  the  system  (3.1)  and  the  relation  (4.6),  we  know  that 

(P(T)h,  h)  -  {Pj (o).  / (o))  +  J  <Dl/f(s),f}T(s)>ds 

O 

+  f  <D2P/f{s),P/?(s)>ds  , 

where  we  have  written  /J7 (0  instead  of  fl  to  emphasize  the  dependence  on  T. 
Recall  that  Pf?  -  -ff. 


t  Of  course  if  <D\(t)v,  v>  2  MvltJ,  (4.7)  holds  as  well  as  the  controllabi¬ 
lity  property  (see  Proposition  3.2).  But  (4.7)  is  weaker. 


'  i 
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In  a  similar  manner  we  can  write  the  relation 

(P(T)h,  h)  «  (/’(/)/(/), /</))  +  ^  <Dlf{s),f{s)>ds 

*  j  (4.9) 

+  J  <D2pf(s),PpT(s)>ds, 

i 

which  holds  for  any  t  €  ( o ,  7).  Therefore  it  follows  that 

ft{P{0p7 (/),/(/))  -  <£>,(*/■  (/),/(/)> 

+  <D2(t)P(t)f? (/),  P(t)f? (t)> 

and  from  (4.7) 

J-(/>(r)/f  (/),/(/))  *  4^(01*  •  (4  10) 

Note  that,  from  Theorem  3.1  one  has 

(P(r)/(r),/(0)  <  p^r(0l2  - 
Therefore  (4.10)  implies 

^•|T1,2(/^T(r)p  2  ^|P,,2(0/(0P  • 

Hence 

(/>y  («),/(*))  <  (P{T)h,  h)'-{k'p)J  <  p  e'{k,p)T  \h?  .  (4.11) 

Therefore,  if 

(^.^)  *  ^tf2. 

we  have 

IflV)!1  s  £  i*i=  . 

From  (4.6)  it  follows  that 

wm  *  i>o  -  e-lk,2p)T .  d 


4.3  Example 

Let  us  turn  to  the  example  considered  in  the  introduction.  We  shall  take 
H  -  L\a) ,  V  -  H\12)  , 

£  -  H  ,  C  -  /  ,  hence  D\  -  /  , 

£  -  LJ(/~)  ,  C  -  y  m  trace  operator . 


.  .  .  . . i»-***K*»-.  .-h/WUifW*  ••  Ursrt»i/.rt/4> 
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The  system  of  optimality  (3.1)  looks  as  follows: 


da 

dt 


aa  +  P  -  o. 


A 

dt 


-  *J) 


-  o. 


6(o)  -  -Pj(o), 


fl(T)  -  A  • 


(4.12) 


The  assumptions  are  satisfied.  Indeed,  the  detectability  condition  is 
satisfied,  since  using  the  stationarity  we  can  apply  Proposition  3.1,  and  we 
have  seen  that  the  system 

o. 

4  */-  *  *• 
r(o)  -  h , 


has  a  solution  in  L?(o,  «>;  VO- 

We  can  check,  at  least  formally,  that  the  operator  P(t )  will  satisfy  the 
equation 


(dz_  dP£ 

•  da, 


)  +  f  P:P;dT 
r 


-  (z,  C)  V  z,  £  in  H'  .  (4  13) 


Representing  P(t)  by  a  kernel 

P(t)h(x)  -  J  P(x,  (.  t)  h(f)  d4 
yields  the  equation  for  the  kernel 

57  -  a Mp  -  a{p  ♦  J  P(x,  Tj,  t)  P(tj,  i,t)  dr  -  <**-£),  (4 14, 

r 

with  the  boundary  conditions 

§£-  (*.  4. 1)  -  o ,  v  *  c  r,  4  in  n . 

(•*.  4>  0  m  0  ,  Vac  12,  (  in  r , 

P(*.  4>  o  -  P(4<  *•  •) . 


(.  0  -  P,(a,  f)  (kernel  of  Pc)  . 
The  observer  has  been  defined  in  (1.5)  (the  writing  is  formal) 


(4  15) 


remark  4.1.  It  is  important  to  notice  that  the  assumption  (4.7)  does  not 
imply  controllability  (the  situation  is  not  similar  to  that  of  ordinary  differ¬ 
ential  equations). 

remark  4.2.  It  is  fair  to  mention  that  (1.5)  leads  to  more  complex  calcu¬ 
lation  than  (1.3),  but  it  is  less  affected  by  disturbances.  Moreover,  it  re¬ 
spects  the  boundary  condition  of  the  initial  system. 
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